Metabolic Engineering Communications

{beta}-galactosidases (BGs) are essential enzymes widely used in the food industry, particularly in the production of lactose-free products. Among them, the BG from Aspergillus oryzae is of industrial relevance due to its activity at acidic pH and moderate thermal tolerance. However, enhancing its catalytic performance remains a key challenge. Traditional enzyme engineering methods are time-consuming and resource-intensive, limiting their scalability. Recent advances in Artificial Intelligence (AI), particularly those based on Natural Language Processing, offer a promising alternative by enabling efficient exploration of protein sequence space and prediction of beneficial mutations. In this study, we introduce an ensemble-based, zero-shot Protein Language Model pipeline that reconciles predictions from six independent models (ESM2 and the five ESM1v variants) combined with a diversity-aware candidate selection strategy. Applied to the BG from A. oryzae, this approach identified beneficial mutations leading to novel enzyme variants with up to a four-fold increase in catalytic efficiency on oNPGal, a two-fold increase on lactose, and, independently, a T338I variant with markedly enhanced thermostability ({approx}80% residual activity after 24 h at 60 {degrees}C), all without requiring supervised fine-tuning on experimental fitness data. Our results demonstrate that consensus across an ensemble of PLMs can efficiently enrich beneficial substitutions in industrially relevant enzymes and substantially reduce the number of wet-lab candidates that need to be screened. Table of Contents graphic O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=106 SRC="FIGDIR/small/726739v1_ufig1.gif" ALT="Figure 1"> View larger version (29K): org.highwire.dtl.DTLVardef@18084f7org.highwire.dtl.DTLVardef@99a102org.highwire.dtl.DTLVardef@19a64forg.highwire.dtl.DTLVardef@1f59cff_HPS_FORMAT_FIGEXP M_FIG C_FIG

Budding yeast Saccharomyces cerevisiae is a workhorse chassis for producing added food and agricultural compounds. However, building multi-enzymatic pathways for these chemicals often requires iterative genomic integration, underscoring the need for efficient, rapid genome-editing tools that can reliably target transcriptionally active chromosomal regions. In this study, to accelerate strain construction, we established a genome-editing toolkit to rapidly engineer eight loci, highly expressed hot-spots, but nonessential genomic sites suitable for stable pathway assembly. Our approach integrates three key design features: (i) selectable markers to enable rapid screening of edited cells, (ii) extended homology arms that leverage the yeast homology-directed repair machinery for robust genomic integration, and (iii) co-delivery of Cas9 and guide RNAs to promote efficient double-stranded DNA breaks at specific integration sites. The sequence independence of FASTOP relies on the release of integration cassettes from integrative vectors, mediated by restriction digestion at two flanking multiple-cutting sites in the integration module to minimize the risk of introducing sequence errors during PCR amplification of the integration cassettes. Following the introduction of a fluorescent reporter cassette, we observed high integration efficiencies across the target sites. We then integrated the biosynthetic pathway of plant-derived flavonoid naringenin into the hot-spots of the yeast genome using the FASTOP toolkit. Our results demonstrated that upon expressing the five essential genes in simple shake flask culture, naringenin production reached 505.7 mg/L, representing a significant (69-fold) increase over previously reported titers for comparable minimal heterologous pathways in S. cerevisiae. Together, the FATSOP toolkit provides a user-friendly platform for reliably modifying hot-spot loci to rapidly construct multi-enzymatic metabolic pathways in S. cerevisiae, while achieving high production levels for high-value food-relevant metabolites.

Fermentation of C1 gases is an emerging technology where waste gases are bio converted into value-added products. This study navigates the gas fermentation potential of Gordonia rubripertincta to produce carotenoids. The crucial carbon monoxide dehydrogenase (CODH) enzyme, necessary for gas uptake by the microbe, was found to be present in G. rubripertincta through blastp on NCBI website. The organism was then used for gas fermentation experiments in a continuous stirred tank reactor (CSTR) in different modes of reactor operation resulting in the production of about 500 mg pigment/g WCW (wet cell weight). Two important reactor parameters, molybdenum content and pH, were optimized for enhanced carotenoid production. Overall, G. rubripertincta was observed to be an efficient candidate organism for C1 gas fermentation. KEY HIGHLIGHTSO_LIGordonia rubripertincta synthesises aerobic carbon monoxide dehydrogenase enzyme. C_LIO_LIIt is a potential gas fermenting microbe that gives carotenoids as product. C_LIO_LIThe gas uptake efficiency of the microbe is more in fed-batch discontinued mode. C_LIO_LIIn FB-D, the resultant carotenoids are 500+9 mg/g wet cell weight (WCW). C_LIO_LIMo/pH of 20 mg/7.0 resulted in highest carotenoids, i.e., 134+41 mg/g WCW. C_LI GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=87 SRC="FIGDIR/small/722808v1_ufig1.gif" ALT="Figure 1"> View larger version (28K): org.highwire.dtl.DTLVardef@8b1185org.highwire.dtl.DTLVardef@2b6f90org.highwire.dtl.DTLVardef@1a9697dorg.highwire.dtl.DTLVardef@14c9dc8_HPS_FORMAT_FIGEXP M_FIG C_FIG

Cyanobacterial natural products are a rich source of bioactive compounds, yet their heterologous production remains challenging. This study investigates the feasibility of expressing the lyngbyatoxin A (LTXA) biosynthetic gene cluster in a fungal host. The lyngbyatoxin biosynthetic genes (ltxA, ltxB, ltxC) were individually cloned and expressed in Aspergillus oryzae NSAR1 under the control of an inducible promoter. Metabolite production was assessed using LC- MS, and transcriptional analysis was performed by RT-PCR. Codon-optimized constructs and precursor feeding experiments were employed to evaluate pathway functionality. No production of LTXA or pathway intermediates was detected upon co-expression of ltxA-C despite confirmed transcription of ltxB and ltxC. RT-PCR analysis revealed truncation of the ltxA transcript, suggesting incompatibility with fungal transcriptional or splicing machinery. In contrast, expression of a codon-optimized ltxC enabled biotransformation of indolactam V to LTXA in A. oryzae, confirming functional expression of the prenyltransferase. These results highlight transcriptional limitations as a key barrier to heterologous expression of cyanobacterial NRPS pathways in fungal hosts, while demonstrating that downstream tailoring enzymes can remain functional. This work provides insights for future engineering of fungal platforms for cyanobacterial natural product biosynthesis.

Odd-chain carboxylates such as valerate and heptanoate are ecologically relevant metabolites and promising platform chemicals, yet the factors leading to their formation during secondary lactate fermentations remain poorly understood. Here, a continuous anaerobic bioreactor was operated for 297 days under mildly acidic conditions to evaluate how lactate:propionate molar ratios shape product spectrum in lactate fermentations. Valerate was the predominant odd-chain product under all conditions, reaching concentrations up to 110 mM, while heptanoate accumulated only at low levels (<10 mM). At low lactate concentrations (10-20 g/L), product selectivity strongly depended on the lactate:propionate ratio. When lactate:propionate ratios were around 1.2 mol/mol, odd-chain products were favored, whereas higher ratios (up to 4.8 mol/mol) shifted metabolism toward caproate and butyrate formation. However, this trend was not maintained at higher lactate concentrations (30-40 g/L; lactate not fully consumed), where odd-chain selectivities remained high even at lactate:propionate ratios of 4.8 mol/mol. Pathway analysis indicated that under high-lactate conditions up to 30% of lactate was redirected toward propionate and acetate formation, likely via the acrylate pathway. Microbial community analysis revealed a stable dominance of Caproiciproducens spp., that could be correlated to valerate production. Overall, this work provides mechanistic insights into the ecology of lactate fermentations and offers a framework for steering product selectivity in engineered anaerobic systems. HighlightsValerate was the dominant product, reaching up to 110 mM. Lactate:propionate ratios drive product selectivities. High lactate concentrations activated in situ propionate formation pathways. Caproiciproducens dominance was associated with sustained valerate production.

Malic acid is a C4 dicarboxylic acid traditionally produced from petroleum and widely used in the food industry. As a sustainable alternative, it can also be produced as a value-added platform chemical from biomass. Previously, the Escherichia coli strain XZ658 was engineered to produce L-malate via the carbon-fixation reductive branch of the TCA cycle. In this study, we further improved this system by relieving allosteric regulation of citrate synthase, addressing redox imbalance, and enhancing malate export. These modifications approximately doubled the L-malate titer in the final strain MO128 compared to XZ658 under simple batch fermentation conditions. The process achieved a high mass yield of 1.2 g malate g-{superscript 1} glucose, highlighting the carbon-fixation capacity of the reductive TCA pathway for fermentative malate production.

Glutathione (GSH) is a tripeptide that plays essential roles in redox regulation and stress responses across organisms. In Escherichia coli, the GSH-specific {gamma}-glutamyl cyclotransferase (ChaC) has been characterized biochemically, yet its physiological role remains unclear. Moreover, ChaC has been annotated as a regulator of the Na/H antiporter ChaA based on its genomic association, although experimental evidence supporting this function is limited. In this study, we investigated whether chaC and its co-transcribed gene, chaB, are involved in sodium transport or GSH metabolism. Gene expression analyses revealed that chaA, chaB, and chaC are upregulated under salt stress. Functional analyses using deletion mutants showed that loss of chaA reduced salt tolerance, whereas deletion of chaB enhanced tolerance and decreased intracellular sodium levels. In contrast, deletion of chaC had no significant effect on salt tolerance or sodium accumulation. Overexpression of cha genes further indicated that chaA, but not chaB or chaC, contributed to salt tolerance. Importantly, overexpression of chaC significantly reduced intracellular GSH levels, whereas chaB overexpression had no effect. These results indicate that ChaC primarily functions in GSH degradation rather than in cation transport, and that ChaB does not participate in GSH metabolism. Our findings clarify the distinct physiological roles of ChaC and ChaB and provide new insight into bacterial physiology regarding GSH metabolism and ion transport in E. coli.

Species of Bacillus bacteria including Bacillus safensis and Bacillus subtilis are finding increasing uses in biotechnology and bioremediation, thanks in part to their metabolic robustness. Malate dehydrogenase (MDH) is at the heart of central metabolism and thus a better understanding of Bacillus MDH proteins could aid in the optimization of these applications. MDH of Bacillus spp. belong to the lactate dehydrogenase (LDH)-like class of MDHs, otherwise known as the MDH3 class. Despite wide prevalence in nature among prokaryotes and archaea, this typically homotetrameric class is understudied compared to the MDH1 and MDH2 classes found in eukaryotes. We therefore recombinantly expressed and purified MDH proteins from two societally relevant Bacillus spp.-B. safensis and B. subtilis-and characterized them biophysically (via Size Exclusion Chromatography-Small Angle X-ray Scattering (SEC-SAXS) and Differential Scanning Fluorimetry (DSF)) and enzymatically (via spectroscopic activity assays). As expected based on their high sequence identity, the two MDH orthologs had similar properties in most regards, including a tetrameric structure and high susceptibility to substrate inhibition. However, we uncovered differences in conditional thermal stability, in addition to subtle differences in enzymatic activity that offer insight into the workings of LDH-like MDH. Summary statementMalate dehydrogenase (MDH) is a fundamental metabolic enzyme, from microbes to mammals, yet comparably little is known about microbial MDH, especially MDH of the tetrameric MDH3 class. We compare the biophysical and enzymatic properties of two such enzymes from the societally relevant bacterial species Bacillus subtilis and Bacillus safensis, offering useful insight with potential biotechnological implications.

Wickerhamomyces ciferrii is a non-model diploid yeast that naturally produces tetraacetyl phytosphingosine (TAPS), a sphingoid base used in cosmetic and dermatological applications. However, its strong preference for non-homologous end joining (NHEJ) over homologous recombination (HR) limits conventional genome editing, while disruption of LIG4, a core NHEJ gene, compromises cellular fitness. Here, we repurposed native NHEJ activity to develop a homology-independent multicopy genome integration platform for W. ciferrii. The platform combines three optimized donor-design features: telomeric end-shielding with two tandem copies of an 11 bp repeat to improve linear donor persistence, a defective URA5 auxotrophic marker to enrich multicopy integrants, and 5'-phosphorylated donor termini to enhance transformant recovery and integration output. These features were consolidated into the platform vector pTdmVU5. As a metabolic engineering demonstration, multicopy integration of LCB1 and LCB2, encoding the two subunits of serine palmitoyltransferase, increased TAPS titer by 2.7-fold. This work converts the native NHEJ bias of W. ciferrii from a barrier to precise genome editing into a practical tool for pathway amplification and establishes a framework for engineering NHEJ-dominant non-model yeasts.

Hydrogenophilus thermoluteolus TH-1 is a thermophilic hydrogen-oxidizing bacterium capable of producing poly(3-hydroxybutyrate) (PHB) from CO2. To redirect carbon flux for producing other useful biomaterials, we disrupted the acetoacetyl-CoA reductase genes (phaB1 and phaB2), which are central to the primary PHB synthesis pathway. Unexpectedly, the resulting {Delta}phaB1B2 mutant still accumulated PHB under autotrophic conditions, reaching approximately 25-35 % of the wild-type level. Furthermore, PHB accumulation in the mutant was significantly restored when fatty acids (butyrate and oleate) were used as carbon sources, whereas acetate and malate resulted in reduced accumulation. These results suggest the existence of a PhaB-independent PHB synthesis pathway. We propose that intermediates from the {beta}-oxidation of fatty acids are converted to (R)-3-hydroxybutyryl-CoA, bypassing the disrupted PhaB enzymes. Additionally, the basal PHB production from non-fatty acid sources implies the involvement of a reverse {beta}-oxidation pathway. This study highlights the metabolic versatility of strain TH-1 for future metabolic engineering.

Bacillus subtilis is an important chassis for biotechnology, but its use in multiplex genome engineering is limited by low natural transformation efficiency. Here, we compared inducible promoter systems for synthetic activation of the competence regulator ComK and evaluated their effects on the comG operon competence reporter and transformation efficiency. Xylose- and mannitol-inducible systems outperformed IPTG-based constructs and shifted 96-99% of cells into a reporter-positive competent state. However, reporter activation alone did not predict transformation potential. Optimization of culture density and induction timing increased transformant yield 45-fold relative to the initial protocol and 2800-fold relative to the conventional Spizizen method. Disruption of native competence regulatory genes did not improve performance and often reduced transformation output, highlighting the importance of endogenous regulatory circuitry. Using the optimized strain and protocol, we achieved co-transformation frequencies of 11-18% and constructed multiplex spore-display libraries containing fluorescent protein fusions integrated at multiple loci. Screening identified strong dual-display combinations and showed that cargo loading depends on anchor protein, integration locus, and genetic background. SscA fusions supported the highest display capacity and promoted synergistic co-display. Together, these results show improvements in natural transformation-based genome engineering in B. subtilis and provide insight into the construction of multifunctional engineered spores.

The facultative methylotroph model organism Methylorubrum extorquens AM1 is a known lanthanide user, which has shed light on the role of rare-earth metals in biochemistry. The characterization of a methanol dehydrogenase (MDH) protein which requires lanthanides as an enzymatic cofactor outlined the question of how these metals are acquired from the environment. It has been proposed that mesophilic organisms as M. extorquens AM1 can produce siderophore-like molecules, which chelate, transport and traffic rare-earth elements into the microbial cell. Therefore, we performed the bioinformatic and chemical investigation of M. extorquens AM1 by using genome mining, the CAS and arsenazo assay, molecular networking and chemical analytical techniques. Our results showed that indeed Methylorubrum extorquens AM1 harbored a gene cluster to produce metal chelators. The chemical analysis confirmed the production of the known hybrid hydroxamate-citrate siderophores schizokinen A and N-deoxyschizokinen A, which are very likely the side products of the transformation of schizokinen and N-deoxyschizokinen. The determination of the lanthanide chelation activity of the schizokinen siderophores series against three different lanthanides (La, Eu and Lu) showed no coordination activity, thus ruling out the involvement of schizokinen siderophores in rare-earth metal transport.

Prodigiosin is a red pigment produced by various bacteria, including Serratia marcescens. Despite its wide and promising range of biological activities, the large-scale production of prodigiosin is currently limited by its high cost and low yields. Here we propose and optimize an innovative, low-cost, peanut-based solid culture medium that enhances the yield of prodigiosin produced by Serratia marcescens. Colorimetric assays revealed that peanut significantly stimulates prodigiosin synthesis. Further HPLC-MS analysis allowed us to unambiguously identify prodigiosin and shows that our medium specifically improves the yield of prodigiosin. Overall, our innovative culture medium could help lower prodigiosin production costs and, ultimately, open new industrial applications.

Streptomyces bacteria produce a variety of secondary metabolites that hold clinical and agricultural value, yet their biosynthetic potential remains unrealized as many biosynthetic gene clusters are not expressed under standard laboratory conditions. Expression of these clusters is tightly regulated, often by cluster situated transcription factors. The TetR family are regulators whose activity is modulated by small molecule elicitors. Although many TetRs have been characterized, elicitors have only been identified for a small fraction of them. This lack of data presents a limitation in our ability to exploit elicitor-regulator pairs for activation of silent clusters and underscores the need for predictive and testable strategies for elicitor identification. In this work, we test the use of sequence similarity networks (SSNs) as a predictor of elicitor identity using the well characterized TetR protein, JadR2, that has a known elicitor, chloramphenicol. We utilized SSNs to identify JadR2 homologs that may also be elicited by chloramphenicol. We developed a heterologous Escherichia coli reporter system in which regulator activity was monitored using an EGFP readout of DNA binding activity. Using this system, we screened JadR2 and four homologs for responsiveness to chloramphenicol. We found that 3 homologs were elicited by chloramphenicol, all of which were formerly uncharacterized. These results demonstrate that TetR-family proteins can share elicitor responsiveness and that SSNs can be used to prioritize regulators for functional screening. This work establishes a genomics-informed and bioinformatics-guided framework for linking elicitors to their regulator, expanding the toolkit for natural product discovery by unlocking regulatory information across Streptomyces.

A challenge in using plant biomass is its highly recalcitrant nature, which makes it economically infeasible to utilize. In natural environments, various microbes, including bacteria and fungi, are reported to decompose plant cell wall materials such as cellulose and hemicellulose, and there may be undescribed microbes that contribute to the degradation of plant biomass. We focused on isolating novel plant biomass-degrading bacteria and screened more than 100 isolates from the Tomakomai experimental forest in Hokkaido, Japan. Among them, one novel Bacillus species was chosen for whole-genome sequencing. Comparative genomics and a carbon source utilization assay indicated that the isolate belongs to a subspecies of Bacillus subtilis, which we named B. sp. TTS1. Glucose, cellobiose, xylose, xylan, mannose, or mannan was used as the sole carbon source in the minimum medium, and the growth of this bacterium was determined. Furthermore, a proteomic analysis of B. sp. TTS1 was performed using culture supernatants from various polysaccharide-containing media. In the present study, several key enzymes involved in plant biomass degradation were identified, namely {beta}-1,4-mannanase and xylanase, and they were highly enriched in all tested polysaccharides.

Rhodoquinone (RQ) is a recently discovered component of the mammalian electron transport chain (ETC) with a high degree of tissue-specificity. Currently, a lack of pure analytical standards limits efforts to precisely quantify its levels using liquid chromatography-tandem mass spectrometry (LC-MS/MS) and interrogate its biochemical functions within mammalian ETC complexes. Here, rhodoquinone-9 (RQ-9) and rhodoquinone-10 (RQ-10), and their isomeric by-products isorhodoquinone-9 (isoRQ-9) and isorhodoquinone-10 (isoRQ-10), were synthesized from ubiquinone-9 and ubiquinone-10 starting materials. Isomers were separated and purified by flash chromatography and structurally confirmed with nuclear magnetic resonance (NMR) spectroscopy. The chromatographic and fragmentation patterns of both the oxidized and reduced forms of these electron carriers were further characterized by LC-MS/MS, establishing signatures for their confident identification in lipidomics studies. LC-MS/MS analysis of murine kidney tissue with RQ-9 analytical standard spike-in corroborate the identity of the endogenous murine RQ-9 and enable absolute quantification of its levels. Thus, we synthesized and purified RQ-9 and RQ-10 analytical standards that will enable absolute quantification in mammalian tissues and in vitro reconstitution studies on RQ-9 and RQ-10 in the mammalian ETC.

Fluorescent reporters cover a wide range of applications in both basic and applied research. Whether a study involves microscopic imaging to study (co)-localization of proteins, FRET, biosensing, or quantifying gene expression, fluorophores are attractive reporter candidates due to their relatively straightforward in vivo readout. For microbiological applications, a wide variety of fluorescent proteins with varying excitation and emission wavelengths, brightness levels, and maturation times are available. Careful consideration is required when selecting from this large suite of proteins, especially when choosing multiple fluorophores. This is further complicated in phototrophic organisms, which exhibit strong autofluorescence, especially towards the red part of the spectrum, effectively eliminating common candidates such as mCherry. In this study, the specific properties and performance of a selection of fluorescent proteins are systematically evaluated against the background of photosynthetic pigment-derived autofluorescence in the cyanobacterium Synechocystis sp. PCC 6803. Specific readouts of different combinations of fluorescent proteins are also analyzed using high-throughput methods, namely plate reader fluorescent scans and single-cell flow cytometry to quantify fluorescence. The ultimate goal is to assess each fluorescent protein with regard to: 1.) Its ability to be discerned from cyanobacterial autofluorescence. 2.) Its compatibility with other fluorophores in this context. 3.) Its overall suitability in cyanobacterial research. Several highly suitable fluorescent proteins for use in cyanobacteria are identified, including mTagBFP2, mNeonGreen and mScarlet-I and suitable combinations, covering nearly the whole spectrum of visible light. This study expands the knowledge and toolset for current and future researchers and uncovers a whole spectrum of possibilities for fluorescent protein selection in cyanobacterial cell biology.

Over a thousand fungal genomes have been sequenced, yet manually curated genome-scale metabolic models (GEMs) are available for only a limited number of species. Moreover, these models have often been developed independently, leading to inconsistencies in namespaces, compartment definitions, and pathway representations that hinder comparative analysis, the systematic reuse of prior curation efforts, and the integration of consolidated metabolic knowledge. Here, we present the Consolidated Fungal Core Metabolism Model (CFCMM), constructed by integrating thirteen published fungal models spanning Ascomycota, Mucoromycota, and both Crabtree-positive and Crabtree-negative yeasts. We harmonized metabolites and reactions into a non-redundant shared ModelSEED ontological space, standardized compartmentalization, and refined gene-protein-reaction (GPR) rules. Using pathway-level visualization and systematic gap detection, we further improved the integrated network through literature-guided curation to correct stoichiometry, stereospecificity, and pathway architecture. Orthologous protein family reconstruction and functional annotation workflows were used to validate and inform GPR associations, with particular emphasis on ambiguous enzyme superfamilies and membrane-associated components. Using the resulting CFCMM, we built high-quality central carbon core models for each fungus and performed flux balance analysis to quantify ATP-yield variation under aerobic and anaerobic conditions, explicitly evaluating scenarios driven by differences in electron transport chain (ETC) composition. Simulations reproduced the expected fermentative yield of approximately 2 mmol ATP per mmol glucose under anaerobic conditions and separated the thirteen fungi into two bioenergetic groups under aerobic respiration based on Complex I status, with predicted yields of approximately 30 versus 22 mmol ATP per mmol glucose. Forcing flux through the alternative oxidase bypass further reduced ATP yields to approximately 12 and 4 mmol ATP per mmol glucose in Complex I-containing and Complex I-lacking fungi, respectively. Collectively, this work provides a manually curated, ModelSEED-consistent, and extensible fungal core metabolic template, deployed in DOE KBase as a resource for automated reconstruction of central carbon core models from any sequenced fungal genome. In addition, the CFCMM provides modular components for developing GEMs with more accurate energy predictions and enables robust comparative analyses of fungal bioenergetics and core metabolic diversity.

Prenylated isoflavonoids are widely distributed specialized metabolites within the Fabaceae and contribute to various characteristic biological activities for both plants and humans. Several aromatic prenyltransferases (PTs) have been identified in Glycyrrhiza species, which are the most widely consumed crude drugs in traditional Chinese medicine. However, these enzymes do not sufficiently explain the structural diversity of prenylated flavonoids produced in the Glycyrrhiza genus. To identify additional novel PTs, we used elicited cultured Glycyrrhiza glabra roots as source material, in which elicitor treatment of cultured roots increased the accumulation of multiple prenylated flavonoids. To identify the responsible enzyme, PT candidates were screened using G. uralensis transcriptomes, currently the sole publicly available transcriptomic resource within the genus, and a homolog designated GgBSPT1 (BSPT; a broad-substrate prenyltransferase) was subsequently isolated from elicited cultured G. glabra roots. GgBSPT1 differed from previously identified Glycyrrhiza PTs in both amino acid sequence and enzymatic properties. GgBSPT1 catalyzed 3'-prenylation of isoliquiritigenin and 6-prenylation of five flavonoids, i.e., this PT displayed broad substrate acceptance across 20 distinct flavonoid structures. Overall, elicited cultured G. glabra roots enabled the identification of a previously unrecognized PT that is functionally distinct from earlier reported Glycyrrhiza PTs. This study provides a new insight into the metabolic plasticity of Glycyrrhiza species and expands the enzymatic toolkit for future metabolic engineering of prenylated phytochemicals by the unusually broad substrate specificity of GgBSPT1. Main conclusionUsing cultured Glycyrrhiza glabra roots, we identified a new prenyltransferase involved in the formation of a variety of flavonoids, thereby revealing novel prenylated isoflavonoid pathways in licorice.

Plants produce specialized metabolites that function in plant defense and as attractants to pollinators and symbionts. One class of specialized metabolites are terpenoids that are synthesized from universal C5 building blocks via activities including terpene synthases, cytochromes P450, and glycosyl transferases. Some terpenes are highly valued for their use as insect repellants, fragrances, antimicrobial compounds, low calorie sweeteners, flavors, and medicines. Low abundance in target tissues, present in complex mixtures, as well as challenging extraction logistics are barriers to economic sustainable production of these compounds from their native species. While heterologous expression of terpenoid biosynthetic genes is feasible, the potential derivation of the products into conjugates via endogenous cytochromes P450 and glycosyl transferases limits this approach. In this project, we used multiplex gene editing technologies to overcome these challenges by creating novel tomato chassis with altered terpenoid biosynthetic capacity in fruit. Excluding central metabolic genes to minimalize impacts on growth and development, we selected 23 known and potential terpene-related genes expressed specifically in the fruit for gene editing. Fruit production and metabolic profiles of three chassis lines with alterations in the major classes of fruit specialized metabolites indicate loss of these genes is tolerated for fruit production. These combinatorial knockouts also showed modulation of native carbon reallocation toward endogenous sinks beneficial for a biosynthetic chassis. Establishing metabolite-modified fruit chassis demonstrates efficient combinatorial editing of entire branches of plant specialized metabolism, facilitating engineering of heterologous terpenes of industrial interest in tomato fruit.